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The color dipole formulation provides an intuitive picture of hard pro- 
cesses in high energy scattering. Most importantly, this approach allows 
one to calculate nuclear effects in a parameter-free way. I review the rela- 
tion between the dipole approach and transverse momentum factorization 
and present numerical results for open heavy flavor production in proton- 
proton and proton-nucleus collisions. The cross section for heavy quark 
production off nuclei is known to reflect gluon shadowing, but is also af- 
fected by higher twist effects and by finite coherence length effects. 

PACS numbers: 24.85. -Hp; 13.85.Ni 



1. Introduction 

Hadroproduction of heavy quarks at RHIC and LHC is of particular in- 
terest because this process allows one to measure the poorly known gluon 
distribution of protons and especially of nuclei ^H]- So far, almost all our 
knowledge of the proton gluon density at low Bjorken-x {xBj) is derived 
from deep inelastic scattering (DIS) experiments at the Desy ep collider 
HERA. However, since DIS measures the quark distribution, information 
about gluons can be obtained only indirectly through the (linear) QCD 
evolution equations. These equations predict a strong rise of the gluon den- 
sity at low X Bj , simply because gluons are vector particles and are therefore 
most likely radiated with small momentum fraction. 

Naturally, the gluon density cannot increase indefinitely because at some 
critical density, gluons will feel each others proximity and recombine 
Since the probability for recombination is proportional to the square of the 
gluon density, this saturation effect leads to non-linear corrections to the 
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QCD evolution equations. Gluon saturation has been the subject of intense 
theoretical and experimental investigations. 

Nuclei have an advantage over protons for the study of saturation, since 
they have more gluons to start with. Therefore, non-linear effects in the evo- 
lution of the nuclear gluon distribution should set in at much lower energy 
than for protons. Moreover, large nuclei have an intrinsic hard scale, the 
so-called saturation scale Q^, which is proportional to the nuclear thickness 
T^(6). This gives rise to the hope that saturation is calculable in perturba- 
tive QCD. 

Non-linear effects are expected when asTA{b)xG{x) ~ Q^, where G{x) 
is the gluon density of the target and is the hardness of the probe that 
measures G{x). (In the case of heavy quark production, ~ mg, where 
niQ is the mass of the heavy quark.) In the rest frame of the nucleus, this 
condition simply means that the mean free path of a probe of transverse size 
l/Q^ is smaller than the length of the medium. Because of this intuitive 
interpretation of saturation as multiple scattering, it is often advantageous 
work in the target rest frame. 

In what follows, I shall employ the color dipole approach to discuss open 
charm and bottom production at RHIC and LHC. The dipole approach 
is formulated in the target rest frame, but all observables are Lorentz- 
invariant, of course. Since the dipole approach is most easily understood 
and most widely used in DIS, I shall discuss the latter process first, before 
turning to the main subject of heavy quark hadroproduction. 



2. Deep Inelastic Scattering 

In the target rest frame, DIS looks like pair creation in the gluon field of 
the target, see Fig. ^ For simplification, I consider only the case of a longi- 
tudinally polarized 7* with virtuality Q^. The transverse polarization does 
not contain any qualitatively new physics. A straightforward calculation of 
the two Feynman diagrams in Fig. ^ yields for the transverse momentum 
distribution of the quark. 



da {'jIp 



T{x,kl), (1) 




where = a (1 — a) + and a is the light-cone momentum fraction 
of the quark. The gg-pair exchanges a gluon with transverse momentum 
Ut with the target. The latter is characterized by its unintegrated gluon 
density J^{x, kj,). Note that Eq. (P) is also valid for heavy flavors. 
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Fig. 1. Perturbative QCD graphs for DIS. The virtual photon (j* ) fluctuates into 
a (virtual) qq -pair far before the target. The interaction with the target, which is 
denoted by gluon (G) exchange, can put this fluctuation on mass shell. 



In the dipole approach, a mixed representation is employed, that treats 
the longitudinal (7*) direction in momentum space, but the two transverse 
directions are described in coordinate {i.e. impact parameter) space. After 
Fourier transforming Eq. one obtains 



da (YlP ^qX) _ 1 



j da j dVidVae'^^^-^^^i-^^^) (2) 



d'^Kl, 2(2^)2 
X '^^,^gg{a,ri)-^tli^qq{a,r2){<yqq{ri,x) +Crqq{r2,x) - (Tqg{fl -f2,x)|. 

Since the fi are conjugate variables to k±, one can interpret fi as the trans- 
verse size of the qq-paii in the amplitude and r2 as the size of the pair in the 
complex conjugate amplitude. These wavefunctions are simply the 7* — > qq 
vertex times the Feynman propagator for the quark line in Fig. ^ and can 
therefore be calculated in perturbation theory. Explicit expressions can be 
found e.g. in [J. 

The flavor independent dipole cross section aqg in Eq. ^ carries all the 
information about the target. It is related to the unintegrated gluon density 
by 

^99-(^ ^) = f / {1 - • 

The color screening factor in the curly brackets in Eq. @ ensures that 
(Tqq{r, x) vanishes oc (modulo logs) at small separations. The dipole cross 
section cannot be calculated from first principles, but has to be determined 
from experimental data. 

In the high energy limit, one can neglect the dependence of the gluon 
momentum fraction x on k,j_ and integrate Eq. ((2j over k± from to 00. 
One then obtains a particularly simple formula for the total cross section, 

X — ^ f f '2 T L ^ 

<7tot i7*P X) = 2_^ da d r ^^;^^^(a,fi) (7qq{r,x). (4) 

T,L •' •' 
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This shows how the dipole formula Eq. Q is related to transverse momen- 
tum factorization, Eq. With a phenomenological parametrization of 
aqq{r,x) (see Ref. [S]), an excellent description of HERA low xsj data is 
achieved. 

It was first realized in I7j that at high energies, color dipoles with a 
well defined transverse separation are the eigenstates of the interaction, i.e. 
the eigenvalues of the elastic amplitude operator are given by the dipole 
cross section. For that reason, it is particularly simple to calculate multiple 
scattering effects in the dipole picture: at asymptotically high energy, the qq 
is formed long before the target, and its transverse size is frozen by Lorentz 
time dilation. In that case, one can simply eikonalize aqq{r), 



-^aqq{r)TA{b) 



(5) 



At realistic energies, however, one has to take the finite lifetime of the 
qq-pair into account 4 . Data on nuclear shadowing in DIS are then well 
reproduced 

Note that even though a parametrization of crqq{r, x) takes into account 
all higher Fock states of the 7*, eikonalization accounts only for the rescat- 
tering of the lowest Fock state, i.e. the gg-pair. Multiple scattering of higher 
Fock components containing gluons is however especially important in the 
case of longitudinal photons and for heavy quarks. These need to be treated 
separately, as will be explained in the following section. 



3. Hadroproduction of Heavy Flavors 

Eq. (jlj is a special case of the general rule that at high center of mass 
energies ^/s, the cross section for any reaction a + N {b,c,.. .}X can 
be expressed as convolution of the light-cone (LC) wavefunction for the 
transition a {6, c, . . .} and the cross section for scattering the color neutral 
{anti— a, b,c. . .}-system off the target nucleon A^. In the case of heavy quark 
production this means that the Feynman graphs in Fig. [21 can be written in 
the form PHUl, 

a{GN ^ {QQ}X) = j\a J [^^^^^(a, r)|' dggcla, r), (6) 

where CTqqC is the cross section for scattering a color neutral quark-antiquark- 
gluon system off a nucleon , 

9 1 

crqqcia, r) = - [crqq{ar) + aqq{{l - a)r)] - -(Tqq{r). (7) 
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Fig. 2. The three lowest order graphs contributing to heavy quark production in the 
dipole approach. These graphs correspond to the gluon-gluon fusion mechanism of 
heavy quark production in the parton model. 

The dipole cross section Ugg is the same as in DIS. Again, the LC wavefunc- 
tions ^ Q^QQ can be calculated perturbatively [10 . 

In momentum space, Eq. @ can be written in A;_|_-factorized form jlUl 
111) . Hence the dipole formulation is closely related to the approach of 
Ref. |2j. However, the dipole approach takes into account only the finite 
transverse momentum of the target gluon. The relation of the dipole ap- 
proach to the conventional parton model is explained in Ref. 

In order to calculate the cross section for heavy quark pair production 
in proton-proton {pp) collisions, Eq. © has to be weighted with the pro- 
jectile gluon density. Numerical results are shown in Fig. |SJ Varying free 
parameters, such as the heavy quark mass mg, leads to the uncertainty 
represented by the space in between the curves. See Ref. P for a complete 
discussion. All data except the recent STAR measurement JHI lie inside the 
uncertainty band. Note that one could increase the upper bound in Fig. El 
(right) by allowing renormalization and factorization scales to be different. 

Multiple scattering effects in nuclei become important at low X2-, when 




Fig. 3. Theory vs. data for open charm production. Left panel: varying free para- 
meters in the dipole approach and in the NLO parton model (NLO PM) JWi gives 
rise to the uncertainties. Right panel: Single inclusive rapidity distribution in the 
dipole approach for different values of the charm quark mass. 
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the lifetime of the QQ-paii becomes large enough to allow for more than 
one interaction with the target. There are two different sources of nuclear 
suppression for heavy quark production _10,. I discuss these now. 

First, the QQ pair in Fig.|21can rescatter several times inside the target. 
Since the amplitude for interaction with a single nucleon is known and 
because dipoles are interaction eigenstates, all these rescatterings can be 
resummed by eikonalizing aqq as in Eq. ©. Note that since the typical 
size of the pair is ~ l/mg, double scattering behaves parametrically like 
Aa/a oc Ql{x2) /niQ. Though formally higher twist, nuclear suppression 
due to heavy quark rescattering is enhanced by a factor of the saturation 
scale in the target nucleus, Ql{x2) oc A^/'^ and cannot be neglected for 
charm quarks. 

In addition, there is the leading twist gluon shadowing ^Ij. One of the 
quarks (or the gluons) in Fig.|21can radiate another gluon, which then prop- 
agates together with the QQ-pair through the nucleus. The magnitude of 
the rescattering correction is determined by the transverse distance re this 
gluon can propagate from its parent quark. Unlike the size of the heavy 
quark pair, re is limited only by nonperturbative QCD effects. In Ref. jl4j . 
these nonperturbative effects were modeled for the case of DIS by an at- 
tractive interaction between the gluon and the quark it is radiated off. In 
order to explain the smallness of the Pomeron-proton cross section, a value 
re < 0.3 fm is needed, resulting in a rather weak gluon shadowing in DIS. I 
argued in Ref. jTH] that the gluon propagation radius of 0.3 fm is determined 
by properties of the QCD vacuum and therefore process independent. 

Gluon shadowing RG{x2,b) at impact parameter b is included in the 
calculation through the modified eikonal approximation, 

a^^(x2,r) = 2 J dHl^l-exp(^-^agg{x2,r)TA{b)RGix2,b)^Y (8) 

Note that the leading twist gluon shadowing has the effect of making the 
target nucleus more dilute, thereby reducing the saturation scale and 
diminishing non-linear effects. 

Results are shown in Fig. |3 The eikonal approximation assumes that 
the lifetime of the heavy quark pair (the coherence length Ic oc I/X2) is 
much larger than the nuclear radius Ra- This is certainly the case at the 
LHC (y^ = 8.8 TeV). For ycc = at RHIC, however, Ic ~ Ra and a more 
sophisticated calculation using the Green's function technique developed in 
Ref. [1] would yield a slightly smaller suppression. Gluon shadowing sets in 
at smaller X2 than quark shadowing, since the QQG-fluctuation has a much 
shorter lifetime than the QQ-pair alone because of the larger invariant mass 
of the QQG-state J^. This is taken into account in this calculation. In 
particular there is no gluon shadowing at ycc = at RHIC. The rescattering 
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Fig. 4. Rapidity dependence of nuclear suppression for open heavy flavor production 
inpAu = 200 GeV) and inpPb (^/s = 8.8 TeV) collisions. Dashed curves show 
gluon shadowing only, while solid curves also include rescattering of the heavy quark 
pair. 



of the heavy quark pair is a sizable contribution to the nuclear suppression of 
open charm production, especially at ^/s = 200 GeV, where gluon shadowing 
is weak. The rescattering of 66-pairs at the LHC is expected to be only a 
small effect. 



4. Summary 

At high energies, hadroproduction of heavy quarks can be formulated in 
terms of the same color dipole cross section as low XBj DIS. Even though 
the dipole formulation is most intuitive for DIS, the cross section for sev- 
eral processes can be expressed in terms of agg{r), e.g. Drell-Yan dilepton 
production 17^ and deeply virtual Compton scattering jl8|| . It is not nec- 
essary that a gg-dipole be present diagrammatically. The dipole approach 
is closely related to transverse momentum factorization, the dipole cross 
section being the Fourier transform of the unintegrated gluon density. 

Using a successful phenomenological parameterization of agg{r,x) [H], 
data on open charm hadroproduction in pp are well described [T], except for 
a recent STAR measurement jl3j . 

The advantage of the dipole formulation is that it allows one to cal- 
culate nuclear effects without introducing additional free parameters. All 
parameters are fitted to scattering off protons. The dipole approach takes 
into account both, rescattering of the QQ-pair as well as leading twist nu- 
clear gluon shadowing. For open charm production, both of these effects 
are important, but for open bottom production gluon shadowing dominates. 
An unambiguous extraction of nuclear gluon shadowing is therefore possible 
from a measurement of open bottom production at the LHC. 
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